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In order to determine the orbital characters on the various Fermi surface pockets of the Fe-based 
superconductors Bao.eKo.4Fe2As2 and FeSeo.45Teo.55, we introduce a method to calculate photoe- 
mission matrix elements. We compare our simulations to experimental data obtained with various 
experimental configurations of beam orientation and light polarization. We show that the photoe- 
mission intensity patterns revealed from angle-resolved photoemission spectroscopy measurements 
of Fermi surface mappings and energy-momentum plots along high-symmetry lines exhibit asymme- 
tries carrying precious information on the nature of the states probed, information that is destroyed 
after the data symmetrization process often performed in the analysis of angle-resolved photoemis- 
sion spectroscopy data. Our simulations are consistent with Fermi surfaces originating mainly from 
the d xy , d xz and d yz orbitals in these materials. 



PACS numbers: 74.25.Jb, 74.70.Xa 

I. INTRODUCTION 

The spectral intensity measured by various experimen- 
tal tools is modulated by matrix elements sensitive to 
the nature of the states probed, as well as to the exper- 
imental setup. For particular configurations, symmetry 
imposes some matrix elements to vanish or to reach max- 
ima. Taking advantage of such selection rules, one can 
extract precious information on the probed states. For 
example, numerous textbooks describe how to use Ra- 
man and infrared selection rules to reveal the symmetry 
of phonons and other excitations. As with other probes, 
symmetry plays an important role in the photoemission 
process. It has been used often in the past to identify the 
nature of the electronic states of various systems [IH7] . 

Unfortunately, simple photoemission selection rules are 
restricted to a few configurations, which are not necessar- 
ily accessible with every experimental setups. Moreover, 
a slight sample misalignment may cause a misinterpre- 
tation of the data. In fact, the intensity variations in 
the momentum space often look strange and asymmet- 
ric, and they are usually neglected by ARPES experi- 
mentalists, which refer to them as the nebulous matrix 
element effects. In some cases, Fermi surface mappings 
are symmetrized to make them look more "natural" . De- 
spite several attempts reported previously to reproduce 
experimental data [7MlO|. the determination of the orbital 
characters in ARPES experiments is still not performed 
routinely, mainly due to the complexity of the calcula- 
tions. A simpler and more practical approach is needed 
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to extract useful information that is otherwise commonly 
sacrificed. 

In this paper, we develop a systematic but simple ap- 
proach to the calculation of photoemission matrix ele- 
ments in Fermi surface mappings. We apply this tech- 
nique to optimally-doped Bao.6Ko.4Fe 2 As 2 , a multi-band 
Fe-based superconductor for which plenty of data is 
available in literature [TT], and to FeTeo.55So.45, an Fe- 
chalcogenide superconductor. A precise knowledge of 
the determination of the orbital characters of the low- 
energy bands is particularly crucial in these materials, 
for which superconducting pairing mechanisms involving 
orbital fluctuations have been proposed |l2j. Our cal- 
culations show remarkable agreement with experimental 
data in multiple experimental configurations of polariza- 
tion and beam orientation. 



II. EXPERIMENT 

In order to test our numerical approach, we performed 
ARPES experiments on high-quality single-crystals of 
Bao.6Ko.4Fe2As2 and FeTeo.55Seo.45 under various con- 
ditions. For each experimental setup, samples have been 
cleaved in situ and maintained in ultra-high vacuum con- 
ditions. ARPES Fermi surface mappings were performed 
at the Institute of Physics, CAS, in a weakly polarized w- 
configuration using a MBS Tl microwave-driven helium 
source {hu = 21.2 eV) and a VG-Scienta R4000 elec- 
tron analyzer. Synchrotron-based experiments were also 
performed at Swiss Light Source beamline SIS and at 
beamline UE112_PGM-2b of BESSY using a VG-Scienta 
R4000 electron analyzer mounted in p and s configura- 
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tions, respectively. For these experiments, photons in 
the 20-138 eV range with different circular and linear 
polarizations were used. All measurements have been 
performed below 20 K. 



III. DEFINITIONS AND CONVENTIONAL 
SELECTION RULES 

Photoemission is a complex quantum problem which 
is far from easy to handle. For a simpler description, 
it is very convenient to decompose this process into the 
three steps of the so-called 3-step model [13J: (i) excita- 
tion of an electron of initial state | i) into a bulk final 
state; (ii) travel of the excited electron towards the sur- 
face; (iii) transmission of the excited electron through the 
surface into a final state | /) approximated by a plane 
wave. During the whole process, the relaxation of the 
remaining electrons and their interactions with the pho- 
toelectron are neglected. Within the 3-step model, the 
matrix element characterizing the photoemission process 
is given by: 



M if = (/ | A • r 



(1) 



where A is the potential vector associated with the in- 
coming photon and r is the position operator. For sake 
of clarity, we also disregarded a constant prefactor. 

We present in Figure [IJa) two commonly used ARPES 
configurations that simplify the analysis significantly. We 
call A^ and A CT , respectively, the components of the po- 
tential vector parallel (it polarization) and perpendicu- 
lar (a polarization) to the emission plane defined by the 
vector k along which the photoemitted electron is ejected 
and the normal to the sample surface. Similarly, the in- 
cident plane is defined by the incident light vector and 
the normal to the sample surface. When dealing with 
unpolarized light, it is also useful to define two special 
configurations of ARPES setup. Hereafter, we call p and 
s the ARPES configurations for which the incident and 
emission planes are parallel and perpendicular, respec- 
tively. With 9i described as in Figure [T| a), the potential 
vector can be expressed in a more general configuration 
with linear polarized light as: 



A = (A x ,Ay,A z ) = (-A^cose h A a ,A„smdi) (2) 

Right-handed circular polarization C + and left-handed 
circular polarization C_ are defined by A(C±) = A n ± 
iA a , and thus for circular polarized light we have: 



A(C±) = (-A Tr cos6 l ,±iA (T ,A K sm6i) 



(3) 



Non-polarized light is treated by adding separately the 
contributions of tt and a linearly polarized photons to 
the photoemission intensity | Mif | 2 . 
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FIG. 1. (Color online) (a) Definitions of the p and s ARPES 
configurations, along with the various angles used in the cal- 
culations, (b) and (c) illustrate, respectively, the even and 



odd combination of the d xz 
definition of the 6fs angle. 



and d yz orbitals, along with the 



Since | Mif | 2 is a scalar observable, it must neces- 
sarily transform under crystal symmetry operations like 
the fully symmetric irreducible representation Ti of the 
corresponding group in order to be different from zero. 
In other words, the decomposition of the tensor product 
of Ti, Tf and r op , which are the representations associ- 
ated to | i), | /) and A • r, respectively, must contain 
T\, which is possible only if their total parity is even. 



The plane wave 



is always an even state 



with respect to the emission plane. With respect to that 
same plane, the operator A • r has an even and a odd 
parity for light polarization parallel (A w ) and perpen- 
dicular (A CT ) to the emission plane, respectively. Know- 
ing the parity of both A • r and the final state from the 
experimental configuration, one can deduce the parity 
of the initial state by choosing a proper set of coordi- 
nates. For a tetragonal system with d electrons like the 
Fe-based superconductors, the most natural orientations 
for ARPES experiments is to align the sample (i) with 
the Fe-Fe bonds parallel to the emission plane to probe 
the electronic states along the T-M direction (here de- 
fined in the 1-Fe/unit cell representation), or (ii) with the 
Fe-Fe bonds at 45° from the emission plane for ARPES 
measurements along the T-X direction. In these simple 
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cases, the five orbital wave functions d z 2, d xz , d yz , d xy 
and d x 2_ y 2 form a convenient basis to describe the initial 
state. It is often preferable though to use linear combi- 
nations of d xz and d yz to construct the wave functions d 
and c? e , which are odd and even with respect to any emis- 
sion plane, respectively, as shown in Figures[TJ(b) and (c). 
More specifically, we have 



where, 



d e = d xz cos 9 F s + d yz sin 9 FS 
d = -d xz sin 9 FS + d yz cos 9 F s 



(4) 
(5) 



where 9 F $ is the Fermi surface angle defined in Figures 
[TJb) and (c). Although such approach has been used al- 
ready to study the Fe-based superconductors [14HT9] , the 
various interpretations are not always consistent, thus 
calling for alternative methods for determining the or- 
bital characters. 



IV. COMPUTATIONAL DETAILS 

In this section, we explain briefly how to use ARPES 
intensity patterns to determine the orbital characters of 
the Fe 3d electronic states near the Fermi level of Fe- 
based superconductors. A more detailed calculation is 
given in Appendix A. Here we focus only on the main 
steps. 

Within the 3-step model, as mentioned previ- 
ously, we use the 3c? atomic orbital wave functions 
{d xy , d xz , d yz , d z , d x i_ y i \ to characterize the initial state 
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m=-2 
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where a m are coefficients, Rz2{r) oc r 2 e~ r / 3 is the lon- 
gitudinal part of the 3d atomic wave functions with r 
given in Bohr radius units and Y [ m (9, <f) is the spherical 
harmonic with angular moment I and azimuthal moment 
m. The final state, here approximated by a plane wave 
function, can be expressed in terms of the spherical har- 
monics as: 



(r | /> = r 
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where ji(kfr) is the Bessel function. The photoemission 
matrix element Mh associated with the different spheri- 
cal harmonic Y 2 X becomes 



M? f oc (/ | At I i;m = A) 

= (A X ?* + A y r$ + A z r x z ) 



(8) 
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(10) 



= j> d9d(j) sin 8Yf 



(ii) 



Px = x/r -- 

Py = y/r = i 
Pz = z/r = Y° 



(12) 



The passage from these matrix elements to matrix el- 
ements involving the 3c? orbital atomic wave functions is 
performed using the following relations: 
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In Figures [2]( a)- (e), we give the cf> dependence of the x, 
y and z components of these matrix elements for a pho- 
ton energy of 21.2 eV, which corresponds to the la line of 
conventional He discharge lamps, and for kf\t = 0.37r/a, 
where a is the in-plane lattice parameter. We used the 
fact that the standard Gaunt coefficients (/, /x) are non- 
vanishing only for I = 1, 3. In addition, we found empiri- 
cally that the coefficients pz(kf) — 1 and pi(kf) = —2/5 
reproduce the experimental data very well over a wide 
range of photon energy. For a better comparison, all the 
matrix element weights are normalized by z — M (d e ) . We 
note that matrix elements for the purely in-plane orbitals 
d xy and d x i_ y i are smaller than the other ones by a factor 
of 5, even though d xy , d xz and d yz are equivalent orbitals 
under symmetry operations. This effect is caused by the 
smallness of the angle 9k when using a photo energy in 
this range. We also point out that the z component of 
the c? z 2 matrix element is larger than any other, which 
indicates that the d z i matrix element is more sensitive 
than others to a A z polarization. 
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FIG. 2. (Color online) (a)-(e) Angular dependence of the x 
(in red), y (in green) and z (in blue) components of the pho- 
toemission matrix elements related to the various 3d orbitals 
(see the text), (f)-(j) Same as (a)-(e) but using the simplified 
matrix elements (see the text). We used hu — 21.2 eV and 

kp = 0.37T. 



Due to the fact that k z is not a good quantum number 
in photoemission experiments, we introduce a few empir- 
ical parameters to the formula describing the full matrix 
element Mf* and improve the agreement between simu- 
lations and experimental data. The full matrix element 
is now expressed as: 



M s if = (A X M S X + A y M 5 y + A Z M 5 Z * w z * e*<*)w s (14) 

where 5 = d z 2 , d xz , d yz , d xy , d x 2_ y 2 . From direct compar- 
ison with experiments, we found out that the ratio be- 
tween the matrix elements associated with the d xy and 
d yz orbitals is only 1/2 instead of 1/5. We thus intro- 
duced the weight factor w$ = {5(6 = d xy , d x 2_ y 2), 1 oth- 
erwise}. Within our semi-quantitative approach, these 
parameters are viewed as phcnomcnological parameters 
compensating for our simplified model. They have been 
fixed at the same values for all our simulations. For the 
study of Fe-based superconductors, w z and 7 have been 



fixed to 4 and k z c+ , respectively, where c is the lattice 
parameter along the z direction. 

Since ARPES allows only measurement of the inten- 
sity lfJ(j>k) — \Mfj\ 2 , and more precisely the relative 
distribution of intensity in the momentum space, several 
prefactors can be dropped in the calculations, including 
the imaginary prefactor. Assuming small 9^ in Eq. (|9j), 
we can simplified the matrix elements as follows, keeping 
only the angular parts of the matrix element components: 

M(d z 2) = (-cos^fc, -sin^-l) 

M(d xz ) = (1,0,-cos^) 

M{d yz ) = (0,l,-sin^ fc ) 

M(d xy ) = (sin0 fc ,cos0fc,-sin2(/>fc) 
M(d x 2^ y 2)= (cos0 fc , - sin^fc, - cos20 fc ) 

M (d e ) = M (d xz ) cos 9 FS + M (d yz ) sin 6 FS 

M(d ) = —M (d xz ) sin 9 FS + M (d yz ) cos 9 FS 

The different components of these simplified matrix 
elements are given in Figures [2[f)-(j). Although their 
precise absolute values differ from that of the components 
in Figures[2|a)-(e), they carry essentially the same orbital 
information while simplifying calculations significantly. 

We now consider the effect of light polarization on the 
photoemission response, which is widely known by ex- 
perimentalists to be important. We first start by the 
experimental observation of a difference, often called cir- 
cular dichroism, between the photoemission response to 
left-handed and right-handed circular polarizations. This 
effect can have different origins 20 . For example, it 
has been attributed to spontaneous breaking of the time- 
reversal symmetry in Bi 2 Sr 2 CaCu208+5 [21]. This effect 
is quite different from the circular polarization used in 
YBa2Cu307_a, for which circular dichroism appears as a 
surface anomaly. It has been useful to separate the pho- 
toemission contributions of the bulk and of the highly 
polar surface resulting from the absence of natural cleav- 
ing plane in this material [2"2~]--24j . In this particular case, 
only a non-trivial combination of the photoemission re- 
sponses to left-handed and right-handed circular polar- 
ized light can allow a full separation of these two compo- 
nents ps] . 

We note that besides these anoumalous circular dichro- 
ism effects, one should also expect asymmetric photoe- 
mission response to C+ and C_ light, depending on the 
geometry of the ARPES configuration. Indeed, it has 
been shown that unless the photoemitted momentum, 
the normal to the sample surface and the incident beam 
momentum are all coplanar (p ARPES configuration) in a 
mirror symmetry plane of the sample, circular dichroism 
can be observed [20]. Rather than searching an origin 
for circular dichroism in the Fe-based superconductors, 
which goes beyond the purpose of the current work, i.e. 
to extract useful information on the orbital characters of 
the bands and FSs observed by ARPES, here we simply 
try to describe its phenomenology and to add it as a tool 
to determine the orbital characters of bands. 
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By working out the details of Eq. (13), one can show 



11,1 i i r^ =d z 2 ,d xz ,dy Z ,d Iy ,d 2 _ 2 

that all the matrix elements Az a = x ,j/,z 
have pure imaginary values. Assuming the form of the 
potential vector given in Eq. ^ for circular polariza- 
tion, we deduce that the photoemission intensity for C± 
polarized light is given by: 



4 ± « 
oc 



\M S C J 2 

\A„(- cos 6iM 5 x + sin 0jM*) ± iA a M y 



■S\2 



(15) 



-sinftM/f| 2 



The previous equation indicates that there should be 
no difference between the photoemission responses to C+ 
and C_ polarized light if all the matrix elements have 
pure imaginary values, which is supposed from Eq. (13 1. 



To account for the difference occurring in real experimen- 
tal data, we add a phase to each matrix element. The 
photoemission intensity when circular polarized light is 
used thus becomes: 



I 5 C cx ^(-cosfl/Mfe* 7 * +sin0;M*e i7 *) ±iA a M*e 



ilv | 2 

(16) 



The addition of phase factors to each matrix element also 
influences the photoemission intensity 1^ corresponding 
to 7r-polarized light and the photoemission intensity I a 
associated to cr-polarized light, which are now respec- 
tively expressed as: 



4 cx \A„(- coBfljATje* 7 - + sin6» z Mfe^)| 2 (17) 



and 



with a previous work [20J . This contrasts with the inten- 
sity predicted for a non-polarized light excitation. Since 
the photoemission intensity I non for a non-polarized light 
excitation can be described by the sum of 1^ and I IJl we 
have: 



'(ttJoc \A a M^y\ 2 

(?r) cx I A % (- cos OiMie** + sinftMfe'' 7 *)! 2 



jS—even 



These equations indicate that in the p ARPES config- 
uration, even symmetry orbitals may lead to an intensity 
asymmetry along k x , but not the odd symmetry ones. 

In the s configuration, electrons are collected in the 
k y — k z plane, and we should expect different selections 
rules. Indeed, we now have for the odd symmetry or- 
bital characters: M 5 x = odd ^ 0, but M 5 z = odd = and 
My =odd = 0. For the even symmetries, the experimental 
configuration imposes M s x = even = 0, but M*= even ^ 
and My =even ^ 0. Consequently, the photoemission in- 
tensity response to circular polarized light in the s con- 
figuration becomes: 



rS—odd 
L C± 



i: 



(a) cx \A W (- cos OiM^)] 2 
l (a) cx |A ir (sm0jM J Je i7 ') ±iA <T M*e i ~ f * 



In contrast to the p configuration, the equations show 
that we can expect circular dichroism in the s configura- 
tion, in agreement with a previous work focused on core 
levels [5DJ. The use of circular polarized light is also a 
useful way to determine the symmetry of the band struc- 
ture. As for the photoemission response to non-polarized 
light in the s ARPES configuration, we now have: 



I s a oc KM*e^| 2 (18) 

These later considerations allow us to predict appro- 
priate phenomenological forms for the photoemission in- 
tensity responses to circular and unpolarized light in 
the common p and s ARPES configurations. As illus- 
trated in Figure [TJa) , the photoemitted electrons in the 
p-type ARPES configuration are collected in the k x — k z 
plane. For odd symmetry orbital characters, we then 
have M s x = odd = 0, M 5 = odd = 0, but M & = odd ^ 0, while 
for even symmetry we get M x =even ^ 0, M s z = even ^ 0, 
but My =even =0. As a consequence, 

4± od< V) « I ± iA a Mye ilv | 2 
I s c =l even {-K) cx \A^{- cos 9iM s x e 1 ^ + shift M 2 V 7 *)| 2 

From these later equations, we can conclude that in 
the p configuration, there is no difference between C+ 
and C_ along the high symmetry line cut, a result valid 
for both odd and even orbital characters and consistent 



It: dd (°) cxK(-cosftMie^)| 2 
In7n ven (*) « K(sinftM z V^)| 2 + \A a M*e^\ 2 

Even though the comparison between the photoemis- 
sion intensity recorded with linear 7r-polarized and cr- 
polarized light give the strongest contrasts, some assump- 
tions on the symmetry of bands can still be made based 
on data recorded with non-polarized light, such as a tra- 
ditional He discharge lamp. 



V. ORBITAL CHARACTERS IN Ba e K 4Fe 2 As 2 

Following LDA band calculations indicating that the 
orbital weight around the Fermi level in Bao.6Ko.4Fe2As2 
is dominated by the Fe 3c? orbitals d xz , d yz and d xy , we 
only considered the related matrix elements in our sim- 
ulations. More specifically, LDA predicts that there are 
three holelikc Fermi surfaces centered at the T point with 
d xy , d e and d Q orbital characters. Previous ARPES re- 
sults also show the existence of 3 holelike Fermi surface 
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pockets centered at T, two of them being nearly degen- 
erate 6, 26-28 . Following a previous notation, here we 
call /3 the outer Fermi surface, and a and a' the two 
others, which will be considered degenerate in our sim- 
ulations. At M = (it, 0), here defined in the 1 Fe/unit 
cell description, theoretical calculations predict a Fermi 
surface pattern formed by the hybridization of 2 ellipses. 
For k z = 0, the ellipse tips have a d xy orbital charac- 
ter while the inner part comes from d yz and d xz [53 |30j ■ 
This orbital distribution around M is reversed for k z = n. 
Theoretical calculations also predict a non-negligible k z 
variation at the M point [2T)H3T] that is not observed 
by ARPES [21]. While ARPES performed for several 
Fe-based materials with different cleaving surfaces reveal 
k z variations of the electronic band structure at the T 
point [llj . the reasons behind this experiment vs theory 
discrepancy for the electronic band structure at the M 
point are still under intense debate. In our simulations, 
we use as M-centered electronlike Fermi surface pockets 
the /c z -invariant hybridized functions determined from a 
three-band model [32] : 

{d X y/yz) =d yz t' 3 ■ ism9 FS 

-d xy (t' 2 -cos 2 FS + e° xy ) (19) 
(d xz / xy ) =d xz t' 3 ■ icos9 F s 

-d xy {t' 2 ■ sin 2 9 FS + e° xy ) (20) 

where we imposed t 3 = t 2 = 1 and e xy = 0.1 for conve- 
nience, these parameters making the weight of different 
orbital characters similar to random phase approxima- 
tion (RPA) results. As we show below, such functions 
are at least consistent with the ARPES observations. 

Figures |3^a)-(d) show the Fermi surface intensity pat- 
terns of Bao.6Ko.4Fe2As2 in various configurations. For 
each experimental pattern, we give in the second column 
from the left the corresponding result from our calcula- 
tions using the orbital configuration given above (Fig- 
ures |3fe)-(h): Simulation A) . The size of each Fermi sur- 
face used in the calculations is chosen to match approxi- 
mately the size of the corresponding experimental Fermi 
surface. We note that small variations in the Fermi sur- 
face sizes do not have qualitative effect on the calculated 
patterns. In the first experimental configuration, light is 
cr-polarized along the x-axis direction. The experimen- 
tal results indicate much stronger weight at the M point 
than for the T-centered Fermi surfaces. The intensity is 
even weaker for the /? band, especially along the k y direc- 
tion. Since the polarization is parallel to k x , this result 
suggests that the /? band must have an odd symmetry 
along both k x and k y , and we thus tentatively associate 
the f3 band to a d xy orbital character, leaving d a and d e 
symmetries for the nearly degenerate a and a' bands. In 
this configuration, our simulation shows a much stronger 
intensity at the M point than for the T-centered Fermi 
surfaces, in agreement with experimental data. More- 
over, it predicts weaker intensity on the outer T-centered 



Fermi surface, with even weaker spectral intensity along 
the k y direction, which is also consistent with the exper- 
iment. 

To test our approach and our orbital attributions fur- 
ther, we show in Figure [3^b) results obtained at 21.2 
cV (near k z = [55]) with light 7r-polarized along 
r-X(7r/2, 7r/2). The Fermi surface mapping is quite 
counter-intuitive, with very strong intensity spots found 
on the T-centered Fermi surface pockets in the first quad- 
rant. The result also indicates strong intensity on the 
tip of the ellipse that has been measured. Surprisingly, 
even such a peculiar Fermi surface pattern is qualitatively 
well reproduced by our simulation displayed in Figure 
(3^f), except perhaps for a weaker intensity on the inner 
T-centered bands than expected. This good agreement 
between simulation and experiment reinforces our initial 
orbital assignment. 

In Figures |3jc) and (d), we present the Fermi surfaces 
obtained with non-polarized light from the la spectral 
line of a Helium discharge lamp (hv = 21.2 eV) for a 
beam incidence aligned along the T-X and T-M direc- 
tions, respectively. Although both configurations give 
rise to much stronger intensity along the M-centered 
Fermi surface elongated along k x than the one elongated 
along k y , the Fermi surface patterns are quite different 
around the Brillouin zone center. While the map ob- 
tained with the T-X orientation of the light shows spec- 
tral intensity almost suppressed in the third quadrant, 
the intensity has a more symmetric distribution in the 
map recorded in the T-M configuration, albeit for an 
intensity slightly smaller below the k x = line than 
above. Moreover, the j3 Fermi surface exhibits an ad- 
ditional suppression of intensity along k x and k y . Once 
more, our simulations, displayed in Figures [3jf) and (g), 
explain well the strange spectral weight intensity distri- 
bution found experimentally. 

In the third column from the left in Figure [3j we illus- 
trate the sensitivity of our approach to distinguish be- 
tween two sets of simulations by displaying simulation 
results (Simulation B) using a wrong orbital assignment. 
As compared to Simulation A, we exchanged the orbital 
characters of the a' (d in Simulation A) and /3 (d xy 
in Simulation A) bands. We also switched the orbital 
characters around the M point, where now the tip is con- 
sidered to have a d xz /d yz character as opposed to a d xy 
orbital character for the inner part. Although the results 
seem also good when using cr-polarized light, the agree- 
ment becomes much worst for other configurations. This 
observation is valid not only for the f3 band, but also for 
the Fermi surface intensity pattern at the M point, which 
is mainly aligned along k y rather than k Xl in contrast to 
the experimental results and to Simulation A. For these 
reasons, we argue that the orbital configuration used in 
Simulation A is at least compatible with the experimen- 
tal results, whereas the one used in Simulation B must 
be discarded. 

Additional information can be obtained from the simu- 
lations away from the Fermi level. In Figure[4]we display 
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FIG. 3. (Color online) Fermi surface intensity patterns of Bao.eKo.4Fe2As2. (a)-(d) Experimental results with different 
photon energies, polarizations and incident beam directions, (e)-(h) Corresponding simulated results using the simplified 
model described in the text (Simulation A: optimized orbital configuration). The inner T-centered a and a' Fermi surface 
pockets with d e and d orbital characters are considered degenerate. The outer one (/3 band) is associated to the d xy orbital. 
The tip of the M-centered Fermi surface pockets has pure d xz or d yz orbital characters while the inner part carries a dominant 
d xy orbital character, (i)-(l) Same as (e)-(h) but using a wrong orbital assignment (Simulation B). The orbital characters 
of the ft and a' bands have been exchanged compared to Simulation A. The orbital characters of the tip and inner part of 
the M-centered FS have also been exchanged. Red double-arrows and blue arrows indicate the in-plane components of the 
orientation of the light polarization and direction, respectively. 



the ARPES intensity plots of Bao.6Ko.4Fe 2 As 2 recorded 
along the M-r direction using 138 eV photons. This pho- 
ton energy corresponds to k z = n, where the a and a' 
bands have the largest separation and thus their appar- 
ent degeneracy is removed [35]. As expected, the inten- 
sity pattern is strongly polarization-dependent. While 
the M-centered bands have very high intensity for a po- 
larization as compared with the T-centered bands, the 
opposite is observed for tt polarization. The spectrum 
obtained with circular polarization is more or less an hy- 
brid of the two others. Interestingly, the spectral weight 
is strongly asymmetric with respect to the zone center 
when using 7r polarization, whereas it is almost symmet- 
ric for the spectrum recorded with a polarization. 

The dispersions and Fermi wave vectors of the various 
bands can be approximated from the intensity plots as 
well as from the corresponding curvature intensity plots 
[33] , which are given in the second row of Figure [IJ Using 
this information, we performed simulations for energies 
away from the Fermi level. The results are compared di- 
rectly to the MDC profiles in Figure [4] To simplify, we 



attributed the same half-width at half-maximum to each 
band. Despite this simplification, the simulations allow 
a good understanding of the MDC profiles. For exam- 
ple, the simulations predict the relative symmetry and 
asymmetry of the photoemission intensity with respect 
to the zone center. More importantly, they allow us to 
pin down the orbital characters of the a and a' bands. 
Symmetry imposes the intensity of the d e band to vanish 
when using a polarization along that particular direc- 
tion. Accordingly, only two bands are observed around 
r in this configuration. In contrast, both the d xy and d 
bands should vanish around T when using 7r polarization. 
Accordingly, only one band is detected around T using 
7T polarization. Since these bands have different Fermi 
wave vectors, their orbital characters appear clearly af- 
ter supcrimposition of the MDC profiles of the spectra 
recorded with the a and tt polarizations, as illustrated 
in Figure |4]j). For instance, we conclude that while the 
innermost band, the a band, has a d Q symmetry, the 
a' band corresponds to the even combination of the d xz 
and d yz orbitals. Our simulations also confirm that the 
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FIG. 4. (Color online) First row: ARPES intensity plots 
of Bao.6Ko.4Fe2As2 recorded along the M-F direction using 
138 eV photons and (a) a, (b) ■k and (c) circular right polar- 
izations. Second row: corresponding ID curvature intensity 
plots [33j along the momentum direction. Third row: corre- 
sponding MDC profile integrated within 10 meV below Ef, 
compared to profiles simulated using the same half-width at 
half- maximum for each band, (j) Direct comparison of the 
MDC profiles recorded with a (red) and n (blue) polariza- 
tions. 



/9 band carries a dominant d xy orbital character. 

Even though the situation is a little more complicated 
around the M point due the weaker photoemission in- 
tensity with 7r-polarized light, our simulations reproduce 
qualitatively well the experimental MDC profiles given 
in Figures |4^g) , Qh) and |4|i) , and suggest that the or- 
bital character at the tip of the electronlike Fermi surface 
pockets with ellipsoidal shape is d xz (d yz ). This conclu- 
sion differs from a previous ARPES study on Co-doped 
BaFe2As2 that rather attributed d x 2_ y 2 and d z 2 charac- 
ters to the tip [TS], which does not show up at the M 
point in LDA band calculations [3TJ [MJ [35] ■ However, 
both ARPES studies indicate that the shape and orbital 
characters of the Fermi surfaces at the M point is pre- 
served along k z , in contrast to LDA band calculations. 

We now investigate circular dichroism for the band 
structure at the T point and demonstrate that it con- 
tains information on the orbital characters of the differ- 
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FIG. 5. (Color online) (a) and (b) give the ARPES intensity 
in the p configuration with C+ and C_ polarizations, respec- 
tively, (c) MDC profiles near Ep corresponding to (a) and 
(b). (d) and (e) Simulated ARPES intensity in the p configu- 
ration with C+ and C_ polarizations, respectively, (f) MDC 
profiles near Ef corresponding to (d) and (e). (g) and (h) 
show the ARPES intensity in the s configuration with C+ 
and C_ polarizations, respectively, (i) MDC profiles near Ep 
corresponding to (g) and (h). (j) and (k) Simulated ARPES 
intensity in the s configuration with C+ and G- polarizations, 
respectively. (1) MDC profiles near Ef corresponding to (j) 
and (k). To simplify the simulations, we chose 7 = 7r/4. 

ent bands. Figures [5ja) and [5jb) show the experimen- 
tal data obtained in the p ARPES configuration using 
C+ and C_ incoming light, respectively. As expected 
from the selection rules derived in the previous section 
for this particular setup and in agreement with our sim- 
ulations displayed in Figures |5jd) and [5^e) , we do not 
observe strong variations between the two sets of data. 
This is also confirmed by the near-i?F MDCs shown in 
Figure p3Fc) as well as the simulated ones given in Fig- 
ure [5jf J7 Interestingly, the experimental data only show 
strong intensity for the degenerated inner band [a (odd 
symmetry) and/or a' (even symmetry)], but not for the 
/3(odd symmetry) band. This behavior is captured by 
our simulations and confirms that the /3 band has a odd 
symmetry orbital character. From our selection rules, 
we deduce that mainly the a' band is observed in this 
configuration. 

The situation becomes quite different for the data 
recorded in the s configuration, once again using C+ and 
C_ incoming light. The corresponding experimental data 
are illustrated in Figures [5](g) and^h), respectively, and 
the MDC profiles near Ep are displayed in[5ji). When 
switching from C+ to C_ polarized light, the observed 
asymmetry in the intensity is qualitatively reversed with 
respect to the T point. As expected from our simula- 
tions given in Figures [5jj) andJSJk), the largest switch 



9 



o.4- (a)i 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 



k x (jt/a) 



FIG. 6. ( a ) _ ( c ) Fermi surface intensity patterns of 

Bao.6Ko.4Fe2As2 recorded with 60 eV circular polarized light. 
The white arrows show the direction of the incoming light. 
(d)-(f) Corresponding simulation results. 
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in the intensity asymmetry is found on the inner band 
that has an even symmetry orbital characters. Although 
observed, this effect is less pronounced for the intensity 
of the (3 band. 

Circular dichroism is also very well illustrated 
by the Fermi surface intensity patterns recorded on 
Bao.6Ko.4Fe 2 As 2 using 60 eV circular polarized light, 
which are displayed in Figures [f|a)-(c). While the in- 
tensity on the right side is much stronger when using 
C+ polarization, the situation is reversed when using CL 
light. This effect is also reproduced by our simulations 
given in Figures [6](d)-(f ) . Interestingly, a comparison of 
Figures [6^b) and (c) indicates that the pattern rotates 
when the beam incidence rotates as well. It is worth 
noting that with circular polarized light the minimum of 
intensity occurs always on one side of the incoming beam 
direction whereas it is observed away from the incoming 
beam side when non-polarized light is used, as suggested 
by Figures [3]^c) and (d). 

At this stage we would like to clarify how we deter- 
mined the phase 75 that appears in Eq. ( 14 ) . Although 



we do not understand its complete meaning, which goes 
beyond the purpose of the current paper, we can intu- 
itively relate this phase to the discontinuity along k z at 
the surface of the sample. To fix this parameter, we mea- 
sured the electronic dispersion along k z , as we now ex- 
plain. Despite the 3D nature of the crystal and elec- 
tronic structures of materials measured in ARPES, this 
technique is so to speak essentially a 2D probe since the 
momentum perpendicular to the surface exposed is not a 
good quantum number. However, within the nearly-free 
electron approximation for the final state [36], access to 
the third dimension of momentum is often possible by 
varying the energy of the incident photons. The momen- 
tum along the z direction is then given by: 



FIG. 7. (Color online) (a)-(e) ARPES intensity plots of 
Bao.6Ko.4Fe2As2 recorded around the V point with C+ polar- 
ized light, (f) Normalized intensity difference between the left 
and right sides of panels (a)-(e). The photon energy has been 
converted into k z using a inner potential Vo = 14.5 eV |28] . 
(g) Photoemission intensity plot of the electronic dispersion 
of Bao.eKo.4Fe2As2 as a function of k z . The data have been 
recorded between 30 and 90 eV using C+ light, (h)-(j) Sim- 
ulations correspondin g to the experimental conditions in (g). 
The phase 7,5 of Eq. ( 14 1 has been fixed to k z c + S. 



k z 



V \hv - $ - Eg) cos 2 6 + V (21) 



where 9 is the angle between the emission direction and 
the normal to the surface, m is the free electron mass 
and Vo is to inner potential, which is determined experi- 
mentally. 

The photoemission intensity is expected to change with 
photon energy due to the photoemission cross section [37] 
and can even show resonances at particular photon en- 
ergies. Photoemission measurements over a wide photon 
energy range can indeed be used to determine the elemen- 
tal characters of the states probes [TJ|s203E]. Experimen- 
tally, additional effects that cannot find a simple expla- 
nation in the photoemission cross section are observed. 
Figures [7]j a)- (e) show such an interesting phenomenon: 
the energy-momentum photoemission intensity measured 
on Bao.6Ko.4Fe2As2 samples with C+ polarized light ex- 
hibits an asymmetry that varies with photon energy. At 
38 eV, the left part of the spectrum has a much weaker 
intensity than the right part. This is no longer the case at 
46 eV, where the two sides show almost equivalent inten- 
sity. The asymmetry is even reversed at 52 eV, with the 
left side of the spectrum being much stronger than the 
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right side. The intensity on both becomes almost equal 
once more at 60 eV before recovering the initial pattern 
at 66 eV. After finding the k z correspondence of each 
photon energy using Vb = 14.5 eV (similar to the value 
reported previously [28] h we can plot the normalized in- 
tensity difference between the left and right sides of the 
spectra as a function of k z . The results are displayed in 
Figure J7jf). Interestingly, the data can be fitted by a co- 
sine function with k z c + 3n/2 as argument, where c = 6.6 
A is the lattice parameter of the primitive unit cell, which 
is equivalent to the distance between Fe layers. 

This strange behavior of the photon energy dependence 
of the intensity extends beyond the 38-66 eV range. Fig- 
ure ^g) reveals oscillations in the 30-90 eV range, as k z 
goes up and crosses different Brillouin zones. This range 
corresponds to k z variations between 7tt/c and ll7r/c. 
The Z positions coincide to k z values with the largest 
kp positions, i. e. k z = 7ir/c, 9tt/c and ll7r/c, whereas 
the r positions coincide with k z = 8tt/c and 107r/c. For 
each Brillouin zone, the signal on the left-hand side is 
quite strong as we increase k z from T to Z, while the 
signal is much weaker on the right-hand side. The sit- 
uation is completely reversed with k z increasing from Z 
to r, with the spectral intensity switched from one side 
to the other. To obtain this effect in the simulations, 
the phase 75 has to be fixed to k z c + ^ over the whole 
range (see Figure Qi)). A variation in the phase leads to 
simulated results completely inconsistent with the exper- 
imental data and justify our choice of phase. However, 
a deeper knowledge of the details of the photoemission 
process would be needed to provide an ab initio value for 
this parameter. 



VI. ORBITAL CHARACTERS IN FeTeo.55Seo.45 

We now check our method to determine the orbital 
characters in FeTeo.55Seo.45. Figures [8]ja) and[8]jb) show 
two Fermi surface intensity patterns of FeTeo.55Seo.45 
recorded with a Helium lamp in the p-type ARPES con- 
figuration. The two measurements differ only by the ori- 
entation of the light momentum, which is aligned along 
T-M and T-X for Figures ^et) and [8jb) , respectively. In 
both cases, the Fermi surface patterns exhibit strong two- 
fold symmetry, with stronger intensity along the light 
momentum direction. According to our simulations, only 
the d Fermi surface can follow this behavior. The corre- 
sponding simulation results for the d Fermi surface are 
given in Figures[8|c) and[8]Jd), respectively. We note that 
although small variations in the Fermi surface size do not 
change qualitatively the simulations, obvious modifica- 
tions appear when the size is modified significantly. For 
example, the large (3 Fermi surface in Bao.6Ko.4Fe 2 As2 
displayed in Figure |3jd) carries the same dominant or- 
bital character as the much smaller /3 Fermi surface in 
FeTeo.55Seo.45 shown in Figure [8|a), which correspond 
to similar experimental conditions. Yet, both the exper- 
imental and theoretical results indicate differences. Nev- 




k x (it/a) k x (at/a) 



FIG. 8. Fermi surface and band structure of FeTeo.55Seo.45. 
(a)-(b) Fermi surface intensity plots recorded with unpolar- 
ized photons (21.2 eV) directed along T-M and F-X, respec- 
tively. (c)-(d) Corresponding simulations for the d band, 
with the light momentum indicated by the white arrows, (e) 
ARPES intensity cut recorded along the blue line in (a), (f) 
Simulations of the intensity cut in (e). 

ertheless, in both cases the Fermi surface patterns exhibit 
a suppression of intensity along the x-axis. 

Figure [8je) shows the ARPES intensity cut along T-M 
for a light momentum aligned along the same direction 
(blue line in Figure [8]ja)). Two bands are clearly ob- 
served, one of them not crossing or barely crossing the 
Fermi level. Actually, a fine study indicates the presence 
of the expected third band, which has a much weaker 
intensity and a kp only slightly larger than that of the 
other band crossing the Fermi level [55] . We display the 
results of our simulations in Figure [8] (f ) for a cut in the 
same configuration, where we assume that the inner band 
carries a d e character while the weak outer one is domi- 
nated by d xy . The main observation is that the d e band 
exhibits a strong asymmetry with respect to T. This is 
indeed what is observed experimentally, reinforcing our 
assumption. We thus conclude that the outer band has 
a d xv orbital character. 



VII. DISCUSSION 

Prior to discuss further the method presented in this 
paper, we would like to comment on the results obtained 
for the orbital characterization of the Fermi surface of 
FeTeo.55Seo.45 and Ba .6K .4Fe2As2. The summary of 
our orbital character attributions for the various elec- 
tronic bands in these materials are displayed in Figure 
[9] For convenience, we spaced the a and a' Fermi sur- 
faces in Ba .6K .4Fe2As2, which are almost degenerate 
in the k z = plane. Except for absolute and relative 
variations of the Fermi surface sizes at the T point, these 
patterns hold for all k z values. We stress once more that 
our experimental observation contrasts with the theoret- 
ical expectation of a switch in the orbital distribution 
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FIG. 9. (Color online) Schematic distribution of the orbital 
characters in Bao.6Ko.4Fe2As2 and FeTeo.55Seo.45- Red: d xz ; 



of the M-centered Fermi surfaces at k z — ir compare to 
k z = [521 130] , which may have important consequences 
for inter-pocket interactions [59] , 

The superconducting gap of Bao.6Ko.4Fe2 AS2 is Fermi- 
surface dependent [26! 23 20] . More precisely, it is about 
12 meV large for all Fermi surface sheets except for the 
6 meV gap found on the j3 band, which carries a domi- 
nant d xy character. The 2A/fcsT" c ratio indicates a pair- 
ing in the weak coupling limit for the /? band. Gaps 
in the weak coupling regimes are also observed for the /? 
band in overdoped Bao.3K .7Fe2As2 [4TJ and underdoped 
Bao.75K .25Fe2As2 [42]. Interestingly, the 2.5 meV gap 
size on the f3 band in FeTeo.55Seo.45 (T c = 14.5 K) leads 
also to a similar ratio [35]. 

From these observations, one could be tempted to ar- 
gue that superconducting pairing is controlled by the or- 
bital character, which for some reason could be less ef- 
ficient for the d xy orbital. However, this argument is in 
contradiction with the observation of a d xy orbital char- 
acter at the M point. In reality, the two electronlike 
ellipses at the M point hybridize and form two distinct 
Fermi surfaces [30] • While the inner one is largely dom- 
inated by a d xy character, the outer one is formed by a 
combination of the d xz and d yz orbitals. Both of them 
show a gap size indicating a strong coupling regime |40j . 
A recent study suggests similar results in FeTeo.55Seo.45 
[38] . Therefore, we conclude that in the Fe-based su- 
perconductors there is no direct correlation between the 
orbital character of a Fermi surface and the gap size. 
Analyses of the gap size on various Fermi surfaces using 
gap functions derived from local antifcrromagnetic ex- 
change interactions rather suggest that the relative size 
of the superconducting gap on a particular Fermi surface 



is determined by its momentum position [28 , 38, 43, 44 . 

The method described in this paper is certainly a reli- 
able and relatively simple way to obtain empirically the 
orbital characters of bands in the iron-based supercon- 
ductors. With a proper choice of basis functions, it can 
be applied to other materials as well. Nevertheless, the 
model has its own limitations. For example, it remains 
quite difficult to determine the orbital characters from 
the Fermi surface patterns in the case of bands with 
mixed characters. Some theoretical assumptions are of- 
ten necessary to guide the analysis. For example, we 
assumed a particular angular distribution for the orbital 
characters of the electronlike Fermi surfaces forming the 
Fermi surface at the M point of Fe-based superconductors 
in order to get a nice agreement between the experimen- 
tal data and our simulations. However, the method is a 
powerful tool to discard some scenarios. 

Another important limitation concerns the determina- 
tion of unknown parameters, such as 75 in Eq. (14 1. 
As explained in Section [Vj we imposed the phase 7,5 by 
looking at the photon energy dependence of the Fermi 
surface pattern. It is clear though that the phase itself 
may carry some important information that is not acces- 
sible directly from our simplify model. From the experi- 
mental point of view, further ARPES studies on different 
materials, involving different electronic orbitals or even 
different transition metals, may help clarifying this issue. 



VIII. SUMMARY 

We introduced a simple method to obtain the orbital 
characters of the various sheets forming the Fermi surface 
of crystals. The method exploits the asymmetries ob- 
tained experimentally in the photoemission intensity pat- 
terns of Fermi surface mappings and energy-momentum 
plots revealed by ARPES in various experimental con- 
ditions of beam orientation and light polarization, in- 
cluding non-polarized light. Our method has been suc- 
cessfully applied to Bao.6Ko.4Fe2As2 and FeTeo.55Seo.45, 
which are two Fe-based superconductors. We showed 
that the multi-sheet Fermi surface of these materials orig- 
inates mainly from Fe 3<i electrons with d xy , d xz and 
d yz orbital characters. Our results suggest that there 
is no direct relationship between the strength of the su- 
perconducting gap on the various Fermi surface sheets of 
these multi-band systems and the orbital characters from 
which they are mainly formed. 
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Appendix A: Details of the matrix element calculations 

In the following, we define: 



(r | /) = e ikf r = to^Mkfr) ]T YTWk, &)*T (9, 4>) (Al) 

1=0 m=-l 

(r\nl o m o )=R nlo (r)Y™°(0,<!>) (A2) 

3 

A • r = (A3) 

a=l 

r a =rd —p a (A4) 

where, a — x,y,z.p x ,p y ,p z are the p-orbital angular distribution functions. Within the 3-step model approximation, 
the photoemission matrix element between an initial state | nl m a ) and a final state | /) becomes 



3 oo |.qo ' {• 

(/ | A • r | nl m ) = ]T A a ]Tz ; • At: / d r]l (k f r)R nla {r)r z ■ £ Y™*(0 fe , fc ) • f d^lY^Y™ (A5) 

a=l 1=0 J° m=-l J 

3 OO I r, 

= E A ° E ll ■ Pi l ° • E Y " 1 * (O^M-f dnYr Pa Y™° (A6) 

a=l 1=0 m=-l J 
3 1 oo I „ 



a=l mi— — 1 1= 



Where we defined p a = X) mi =_i a mi Y" 11 . 

We then use Wigner's formalism for the 3j-symbols: 



/ dnY, m Y^Y, m ° = \ + 1} ' 3 ' + ^ (n n n)( ? 1 ^ ) = W( 1 1 S (A8) 

J i i t y v?i mi mo m mi mo 

which may be non-zero only if Iq ^ 0, m — —(mi + mo) and I = Iq — 1 or Iq + 1. 
Hence, 

<( dnY, m Y" li Y, m " = s?r { r +mi) w( l 1 ^° ) + c+r o+mi) ^( z 1 S (A9) 

J 1 1 <° »>'o-i v m mi m i.'o+i ^ m uii mo 
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Substituting the previous result into equation (A7), we get: 

3 l 



(/ | A • r | nl m ) = ^ A a a ™i ( Al °) 

i—ot nil— — 1 

• [« ,o - i * 1 (*/)(-ir^irx^^,w^(_ ( ^- 1 roi) mi t ) ( An ) 



+ ^ + vSi(fc/)(-i) mo+mi c + i +mi (^'^)^ (- { t c + 1 ' '" 



-(toq + mi) wi toq 



M^ 2 (k) 


= M«(2,Q,k) 






Af*»" (k) 


= i^(M a (2, 


-l,k) 


+ M a (2,l,k)) 


M Q rf -(k) 


= y|(M Q (2,- 


-l,k)- 


- M (2,l,k)) 


Af*»»(k) 


= isj\{M a {2, 


-2,k) 


-M a (2,2,k)) 




= J\{M a {2,- 


-2,k)- 


hM a (2,2,k)) 



(A12) 



^ a mi (j!(mi,i ,rao,k) + (toi, Z , to , k)) (A13) 

a— 1 mi — — 1 
3 1 

^-4 Q a mi G„(TOi,/ ,TOo,k) (A14) 

a— 1 mi=-l 
3 

^A a -M n (a,Z ,m ,k) (A15) 



Where, k= (kf,9 k ,(f) k ) and 

Ma;(lo,mo,k) = M„(l,/ ,TO ,k) = Y^(G„(-l,Z ,TOo,k) - G n (l,i ,"io,k)) 

M y (7 ,TO ,k) = M„(2,/ ,TO ,k) = iJ^(G n (-l,lo,m ,k) + G„(l, l , mo, k)) 

M z (7o,TO ,k) = M„(3,Z ,TOo,k) = G„(0, l , ™o, k) (A16) 

We now apply above results to the case of 3d electrons, for which l — 2. We obtain: 



M a x -» (k) = y-(M Q (2,-2,k) + M a (2,2,k)) (A17) 
The even (d e ) and odd (d Q ) combinations of these matrices for the d xz and d yz orbitals around the T point are 

M*° (k) = cos (k) + sin faMfr- (k) 

Mt (k) = - sin 4> k Mt* (k) + cos 4> k M*>» (k) (A18) 
We can express the previous results by defining: 

y, m (0, <f)) = Cl n Pr{6)e im < 1 ' (A19) 
where P '[""(0) contains all the 9 dependence and GJ" contains all the numerical prefactors. For example, we have: 

^i(0, <j>) = ^ J— • sin 6»(5 cos 2 9 - 1) • 

8 V 7T 

= G^-P 3 1 (0)-e^ (A20) 
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For 3d electrons, this leads to equations (A21)-(A25): 



M x (d z 2) = -i\f2C\Cl 



M y (d z 
M z (d z 



-^ Pl {k f )P^ + \ Pz {k f )Pl 



{-COS (f> k } 

{-sin cj> k } 
{-1} 



(A21) 



d x 



M x {d xz ) 
M y (d xz ) 
M z (d xz ) 



iC\C\ 



2 -pi{k s )P° + gP 3 (fc/)P 3 - p 3 (k f )P! • cos20 fe 



{1} 



- iC\C\ [-p 3 (k f )Pi] sin20 fc ~ {0} 



-\pi{k S )Pl + \ps{k f )Pl 



{-COS (j) k } 



(A22) 



M x (d yz ) = iC\C\ [-p a (*/)if] sin20 fe ~ {0} 
M y (d yz ) = iC\C\ 



\px{k s )Pl + ^ P3 (k f )P° + p 3 (k f )Pi ■ cos2<fe 



{1} 



M z (d 



■\pi(k f )Pl + \p 3 (k f )Pl 



{-sin cj) k } 



(A23) 



i.rij 



M x (d X y) 



-iC\C 2 
-iC\C\ 



\px{k s )P\ + \ P3 {k f )Pl - p 3 (k f )Pl(4 
\pi{k f )Pl + \ P3 {k f )Pl + P3 (k f )Pl(A 



My{d X y) 

M z {d xy ) = iV2C°Cl [p 3 {k f )Pi] {-sin 20J 



cos 2 (j) k 



COS 2 (j)k 



1) 

3) 



{sin (f> k } 
{cos <j> k ] 



(A24) 



M x {d, 



My(d x 2_y2 

M z (d x 2_ y 2 



\pi{kf)Pl + \p 3 {k f )Pl - p 3 {k f )Pl(Aco S 2 ^ k - 3) 
\pi{k f )Pl + \ P3 {k f )Pl + p^P^cos 2 ^ - 1) 



—iC\C 2 
—iC\C 2 

iV2ClCl [ P3 (k f )Pi] {-cos20 fc } 



{cos(f> k } 
{-sin (j> k } 



(A25) 



The main idea behind our simplified approach is to neglect the 4>u dependence of the prefactors preceding the curly 
brackets in the previous equations. We note that M x (d yz ) and M y (d xz ) are set to zero since the terms preceding the 
sin 2<pk function are vanishingly small compared to the other matrix components. The curly bracket terms correspond 



to the components of the simplified matrices given in Section IV 
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